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Abstract 

We consider the gravitational field of infinite straight and stationary 
twisted Baby Skyrmion cosmic string. Using the approximate solution of 
Einstein equations, it is shown that the internal phase rotation (twist) 
along the string axis is responsible for a long-range gravitational accel- 
eration resembling that of massive cylindrical shell. We also study the 
stability and gravitational field of circular loops. When the loop radius 
becomes comparable with the string width, the rigidity energy tends to 
stabilize small loops against radial collapse. The nucleon scale-toroidal 
knot with Hopf charge Q = 1 is found to decay very rapidly on the scale 
of the age of the universe due to low energy cost to fiux lines crossings. 
Such knot is therefore excluded from the dark matter scenario of Spergel 
and Steinhardt. However, the Q = loop, stabilized by rigidity, could 
be a candidate for this scenario. In contrast, the electroweak strings are 
prevented from intercommuting due to much larger energy cost to inter- 
section. This makes them a possible candidate for the solid dark matter 
scenario of Bucher and Spergel. 

PACS number(s): 11.27.4-d, 12.39.Dc, 98.80.Cq, 3.75.Lm 



1 Introduction 



A closed twisted Baby Skyrmion is a topological soliton in a 0(3) sigma model 
including the fourth order derivative Skyrme term. That a stable soliton of this 
kind may exist was suggested some time ago by Faddeev [1]. This model is 
widely known as a Skyrme-Faddeev model. More recently, it has been proposed 
that this model arises as a dual description of strongly coupled SU{2) Yang- 
Mills theory [2], in which the knotted strings may represent glueballs. The 
gravitational field of such glueballs may be of interest in view of the possibility 
that they can be a source of dark matter [3] . This is also suggested by previous 
investigations of cosmological implications of vortons [4]. These are loops of 
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superconducting cosmic strings [5] stabilized against collapse by angular mo- 
mentum of its charge carriers. We note that cosmic vortons in electroweak scale 
were suggested as a source cold dark matter in Ref . [6] . 

In the present paper wc focus on the long-range gravitational force which a 
twisted Baby Skyrmion string exerts on surrounding matter. We first calculate 
the active gravitational mass of a straight string (Sees. 1-6). The circular 
loop formed by bending the straight string is studied in Sees. 7 and 8. The 
gravitational field of the loops is studied in Sec. 9 and Sec. 10 discusses possible 
role of these loops in the cold dark matter problem. 

Recently [7], wc have studied the effect of internal time-dependent phase 
rotation on the gravitational properties of a nonlinear sigma-model with the 
single-axis anisotropic potential. The precession induced by this potential pre- 
vents the lump from collapsing to a point. This method of stabilizing the lump 
has been proposed by Leese [8] in a work yielding analytical solution for the 
fields in the CP^ formulation. Evaluating the energy-momentum tensor with 
the use of this solution, and solving the Einstein equations in the weak field 
approximation, we find that the time-dependent internal space rotation is re- 
sponsible for a long-range gravitational acceleration similar to that of a rotating 
cylindrical shell. Vcrbin and Larsen [9] studied previously this problem numer- 
ically beyond the weak coupling approximation and found that, owing to the 
finite frequency of spinning, the metric 500 is not asymptotically flat consistent 
with the long-range gravitational potential found analytically in our paper [7] . 
For a string, situated parallel to the a:"^-axis, this result is expected by recalling 
the breakdown of Lorentz-boost invariance in the x^— direction caused by the 
time-dependence of spinning-string fields [7,10]. 

Obviously, this invariance can be also lost if the string fields acquire, while re- 
maining time-independent, an a;^— dependent phase factor. The resulting metric 
goo is then expected to yield a long-range gravitational potential. The numerical 
results of Verbin and Larsen [9] are consistent with this expectation. 

Our string model is obtained by trivially extending the 2+1 dimensional 
Baby Skyrme model [11] to 3-1-1 dimensions. The static version of this string 
model has been thoroughly studied in connection with the two-component Bose- 
Einstein condensates by Cho et al [12]. In the time-independent version of 
the two-dimensional Baby Skyrmion model, the size instability of the lump is 
prevented by the presence in the Lagrangian of the Skyrme term which is fourth 
order in the field variables. For the role, this term plays in Derrick's scaling 
argument for stability, we refer to the book by Manton and Sutcliffe [13]. 

2 Lagrangian and Stress-Energy Tensor 

We write the Lagrangian density in the form [12] h = c = 1 



£ 



where fx, v 



0,1,2,3. The basic field is a triplet (j) = (^^ij 02,03) of real 
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scalar fields satisfying the constraint 0.0 = 1. The first term in Eq. (1) is 
the Lagrangian density of the pure sigma-model string. The second term is the 
planar analog of the three dimensional Skyrme term. We note that, in contrast 
with the original Baby Skyrmion model [11,13], the anisotropic mass term is 
not included in Eq. (1). As shown below, in the present model a mass term is 
generated by the twist of the field. For a straight string along the x^— axis, we 
make the following Ansatz [12] for the field 

sin /(r) cos(ni^ + nikx^) \ 

sin/(r) sin(niyj + mfcx"^) I (2) 
cos f ir) J 

where m and n are integers, 27r/fc is the period in the a;^-direction and 
r = [(a;i)2 + (a;2)2]i 

The stress-energy tensor is obtained from the formula [7] 



d 



iV\9\£) 



(3) 



Using the Lagrangian density (1), the (right hand side) rhs of Eq. (3) yields 



(4) 



The metric respecting the cylindrical symmetry of the string, and the twist 
along the a;^-axis is written as 



+2 



g3i{x)dx^dx^ +g32{x)dx^dx^) + g33{r){dx^f 



(5) 



The components of the metric tensor are obtained by solving the Einstein 
equation [14] 



Rfiv = —SttGITh^ — -g^uT 



(6) 



where is the Ricci tensor, T^j, is the stress-energy tensor and T = T^. 

In the weak field approximation, the stress-energy tensor is obtained from 
Eq. (4) by using the fiat space metric r/^^ = diag (1,-1, —1, —1). Defining the 
quantity 



(7) 



we have 



too 



-£ 



Til = 2^'[(5i<^)' + «'(^12 + ^13)] + £ 



(8) 
(9) 
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T22 = lv^[{d2$f + K^{S21 + S23)] + £ (10) 
1^33 = lv^[{d3$f + K^{S3i + 532)] + £ (11) 

3 Active Gravitational Mass 

We first consider the 00-component of Eq. (6). We write = rj^i, + /i^y. 
In tlie weak field limit, /i^jy << ^j^j,, we have 

^^00 = -\{df + di)hoo = - ^ V' hoo (12) 

On the rhs of Eq. (6), we need an expression for r = Too — ^T. Using Eqs. 
(1) and (8-11), we obtain 



^2^2 



p(Si2 + 513 + 523) (13) 



4 

With the Ansatz (2), we obtain from Eqs. (7) and (13) 



r{r) = lv''^'{^ + m'k'){f'smfr (14) 

where /' = df /dr. Wc note that this result is correct only in the presence 
of the twist. If mk = 0, the model described by Eq. (1) is unstable unless it is 
augmented by an anisotropic mass term (see Ref. [13]). 

Using Eqs. (6), (8) and (14), the Einstein equation for hoo becomes 

V' hoo = 167rGT(r) (15) 

This equation can be cast in the form of a Newton-Poisson equation for the 

gravitational acceleration g = —\y(f)g, where the gravitational potential (j)g is 
given (pg = ^hoo [7]. Thus, we get from Eq. (15) 

V.g{r) = -47rGp„(r) (16) 

where Pa{r) = 2T(r) is the active gravitational mass density. Equation (16) 
can be solved for the radial component g{r) of g{r) using the Gauss' law yielding 

g[r) = [ T{r)rdr (17) 



The profile fimction /(r) satisfies a nonlinear differential equation which is 
obtained by minimizing the static energy 



e = 2TT Too{r)rdr (18) 
Using Eqs. (1), (2) and (8), we have 
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roo(r) = \rf [/'^ + + m'k') (sin^ /) (l + K^f')] (19) 
This leads to the following variational equation for the profile function 



2/" 







(20) 



In what follows, we are interested in the asymptotic behaviour of the function 
(17) for large r. If the function T(r) exhibits an exponential decay, then for 
r » Ijj where Id is the decay length, the upper limit of the integral on the rhs 
of Eq. (17) can be replaced by infinity. In this case, the acceleration g{r) ~ 1/r 
which corresponds to the long-range gravitational field of a massive rod. Owing 
to the boundary condition /(oo) = 0, the values of /(r) for large r are small 
and satisfy the linearized form of Eq. (20) given by 

+ rf - {r? + n?}?r^)f = (21) 
The solution of this equation is the modified Bessel function of order n 

/(r) ~ Kn{mkr) (22) 
The leading term of the asymptotic expansion of this function for large r is 

A 



f{r) ~ -j= exp(— mfcr) 



(23) 



Hence, the decay length Id ~ 1/mk and the asymptotic form of Eq. (17), 
valid for r >> 1/mk, is 



9{r) 



2Gma 



(24) 



where 



THa ^ 4:TT J T{r)rdr ^ TTnf j y— -\-w?k'^j [J ' sin fj rdr (25) 

This expression for the active gravitational mass can be rewritten in a form 
showing that nia is proportional to m^k^. In the next section, we establish the 
following condition for the size stability of the string 

(^f2 gj^2 jy^^ ^ ^2^2 j^gi^2 ^y^^ (26) 

Using this relation in Eq. (25), we obtain 

ma ^ nrfm^k^ (l + K^f^^ r sin^ fdr (27) 
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4 Spatial Rescaling and Stability Condition 

Starting from Eqs. (1) and (8), we write the energy (per unit length in the 
direction) 



Tood^x 



d'x 



(28) 



We now consider the transformation of this expression under a spatial rescal- 
ing in the x^ , a;^-space [13] 



a:^ — )■ j/^ = Lix^,x^ — > = ux^,x^ — ?• j/^ = x^, d^x = -^d^y 
This implies a following transformation for the fields 



dy^ 



For the Skyrme terms, we obtain 



dy^ dy"^ 



9y3 



and 



'5i3^/z25i3(y) 



(29) 



(30) 



(31) 



\S2i^li^S23{y) 

Applying the maps (29)- (32) in Eq. (28), we obtain the rescaled energy 



(32) 



e(M) = ^ 



(33) 



From the condition that energy is stationary under spatial rescaling, we 
obtain 



&(m) 
which implies 



= y {ds$fd^y + 2/.«2 j s^2{y)dh = (34) 



j (d3$yd^y = K^ J (di$xd2$yd^y 



(35) 



Using Ansatz (2) for the fields, the stability condition (35) takes the form 
given in Eq. (26). 
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5 Relation of String Width to Twist Rate 



A qualitative understanding of the profile funcion /(r) can be obtained by con- 
sidering Eq. (20) for r — >■ 0. Introducing the Ansatz 

/(r) = TT - CrP (36) 
into Eq. (20) and keeping only the leading terms for r ~ 0, we have 

C(n.^ -p'y"-^+Cm^erP - 2>,^C'n^ (p^ ~ py^'P"* = (37) 

For n = 1, we obtain p = n = 1. In what follows, we confine ourselves to 
the case n = 2, since it allows us to make further progress without resorting to 
numerical work. In this case, the first term in Eq. (37) implies p = n = 2. The 
vanishing of the sum of the second and third terms yields for p = 2 

= S (38) 
Hence, the profile function for n = 2 goes as 

f^r)=n-^r' (39) 
This suggests that the width A is related to the twist rate mk as follows 

To substantiate this prediction, we make a variational Ansatz for /(r) of the 
form 

/(O = cos-^ (41) 

where A is the variational parameter to be determined from the stability 
condition (26). We note that Eq. (41) describes the n = 2 lump solution of 
radius A in the sigma model. Expanding the function (41) into Taylor series 
about r = 0, we have 

2r2 

f{r)^^-— (42) 

Comparing this result with the Eq. (39), we obtain A^ = %K,/mk. 

We now turn to the variational estimate of A. According to the stability 
condition (26), the soliton has a preferred size at which the contribution of the 
Skyrme term is equal to the effective anisotropy energy induced by the twist. 
Introducing the Ansatz (41) into Eq. (26), and performing the r— integrations, 
we obtain the following results 

« " / ^2 ^ dr = ^^ (43) 
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and 



mk rsin^/(ir= (44) 

Jo 2 

We see that the Skyrme contribution (43) goes as A~^. Thus it favors large 
values of A. On the other hand, the anisotropy term (44) favors small A. Eq. 
(26) describes the balance between these tendencies yielding a relation (valid 
for n = 2) 

A^ = ^ (45) 
mk 



For n = 1, the profile function is 



„2 



/W = co«"';:^^ (46) 

For this function, the effective anisotropy energy diverges owing to the power 
law decay for large A. 

Let us estimate the expression (27) for n = 2. Using the profile function 
(41), we obtain 

rua = 'K'Tfrrfk^ + j (47) 

Now, A is not a free parameter, rather it is fixed to the preferred value given 
by Eq. (45). Using this value in Eq. (47), the final expression for the active 
gravitational mass per unit strength of the string becomes 

rUa = W \T^V2\u\ + -u^j (48) 
where we introduced the dimensionless parameter 

u = Kmk (49) 



6 Limit of Pressureless String 

From the expressions (48) and (49) it is apparent that a nonzero value of nia is 
obtained only if the twist rate mk and the magnitude of the Skyrme term k are 
both nonzero. We now show that the magnitude of string pressure decreases 
on increasing the parameter u. For a certain critical value of u, the pressure 
vanishes and the string turns into a pressureless gas in analogy with the vortons 
[15]. 

For a string situated parallel to the x^-axis, the pressure in the a;^-direction 
is given by 

'nsd^x (50) 



8 



Using Eqs. (1) and (11), we have 



^2 r 



^33 = 



{dl(l>f - {d2cl>f + K^Sls + S23) - k'5i2 + {d3cPf 



(51) 



Performing the transverse integral of this equation, and using the size sta- 
bihty condition (35), we obtain 



(52 



(52) 



Since the expression k^(5i3 + '5'23) is proportional to u^, we see from Eq. 
(52) that \p\ decreases as u increases. 
Using the Ansatz (2), we have 

J [(di$y + (d2$y](fx = 2t: + ^sinV)r-cir (53) 



and 



(.S13 + 5*23) d^a; = 277^2 y (/'^ sin V) r-dr (54) 

Introducing Eqs. (53) and (54) into (52), we obtain the condition for van- 
ishing pressure in a general form 



9 • 2 /• 
n sm / 



^ rdr = j (^f ^ sin^ rdr 



(55) 



Confining ourselves to the n = 2 case, we use in Eq. (55) the expression (41) 
and obtain the critical value of u for which the pressure vanishes 



V3 



(56) 



Introducing this quantity into Eq. (48), we obtain the upper bound for the 
active gravitational mass 



niaiu = Ucrit) = 49?7^ 



(57) 



It is interesting to compare the expression (48) and (57) with the M-dependence 
of the energy, per unit length of the string, e, given by Eq. (18). Using the size- 
stability condition (26), we obtain from Eqs. (18) and (19) 



e=^- [ + ^ sin2 / + 2m^k^ sin^ / + K^m'^ef^ sin^ f)rdr (58) 

^ Jo ^ f ' 

Performing the transverse integrations with the n = 2 profile function /(r) 
(given in Eq. (41)), we obtain 
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e = rf {a-k + V2-n'^u + y"^) (59) 

For u — Ucrit — VS, this equation yields the critical energy €criL correspond- 
ing to the case of vanishing string pressure. It is interesting that for this value 
of u, the critical active gravitational mass nia is, according to Eq. (57), exactly 
equal to Ccrit- This conclusion follows more directly by considering the sum 

2 2 f 

^+P=^J (^12 + + S23)d''x = ma (60) 

whore the first equality is the sum of expressions (28) and (52), and the 
second is the consequence of Eqs. (13) and (25). 



7 Stability and Radius of a Circular Loop 

Consider a circular loop of radius R much larger than the thickness of the string 
A. The loop is created by bending the straight string and connecting smoothly 
the periodic ends. A large loop is stabilized by the Skyrme energy and the 
twist along the x'^-axis. This leads to a repulsive centrifugal stress opposing the 
centripetal force due to the tension. For small loops (i?/A ~ 1), the ridigity due 
to elastic bending energy provides an additional loop stabilization mechanism. 
Discussion of this effect is postponed to Sec. 8. 

To study the classical stability of the loop, we confine ourselves to the n = 2 
case and use Eq. (59) to write the energy of the loop in the form 

r /~ 47r / \ ^1 

eioop = 27ri?e = 27ri?r?^ Atv + -n'^\/2mkK + — ymknj (61) 

where the definition u = mkn is invoked. To proceed, we need to relate the 
linear momentum mk to the radius R of the loop. From Eq. (2), we see that the 
transverse fields pick up a phase factor mk per unit length along the coordinate 
x^. Thus, if we denote by •& the phase factor for length x^, we have mk = d^i}. 
Next, we relate the winding number density dsi^ to the topological number A'' 
defined as [5] 

AT = i- y dial's (62) 

For a circular loop of radius R, we obtain from Eq. (62) N = Rdsd yielding 

N 

mk= — (63) 
R 

Introducing this relation into the Eq. (61), we obtain the loop energy as a 
function of the loop radius 

eioop = SttW (R + + K (64) 
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where K is the i?-independent contribution given K = The 
expression (67) has a minimum at the preferred radius Re given by 

Rl = ^ (65) 
Using Eqs. (63) and (65), the equilibrium value of the parameter u becomes 

We = ^ = \/3 (66) 

He 

Recalhng Eq. (56), wc sec that Ug = Uerit- Thus, the equihbrium value of 
u that ensures the stabihty of the ring is equal to the critical value u for which 
the pressure p vanishes. 

Wc now consider stability of a loop of small radius Re . In this case, there is 
a problem that is revealed by considering the ratio Re/X. Using Eqs. (45) and 
(65), we obtain Re/X — 0A5N. This indicates that, for N = 1, the loop tends 
to collapse under the centripetal force due to string tension. 



8 Rigidity of Small Loops 

In what follows, we show that an enhancement of Re/X is obtained from a 
model in which the string energy is augmented by a contribution due to rigidity. 
To estimate this contribution, we follow the analogy with elasticity theory of 
bending rods [16]. Thus, we pictm-e the rigid finite width-string as a bundle of 
filamentary lines each representing a flexible zero width-string. 

Due to bending, lines on the convex side of the string are extended, whereas 
those on the concave are compressed. These deformations are described by the 
relative extension {dx^ — dx^)/dx^ = U33 where U33 is the strain tensor. For a 
simple extension, we can use Hooke's law to find the stress tensor[16] 

(^sair, (f>) = 7oo(r)M33(r, </>) (67) 

where, r, </> are the polar coordinates in the cross section of the string. The 
elastic energy per unit volume of the string is 

V{r, 4>) = = \Tm{r)ul^{r, ct>) (68) 

Integrating this expression over the cross section, we obtain the bending 
energy per unit length of the string 

£6 = -/ rdr dcjiToo{r)ul3{r,cj>) (69) 
^ Jo Jo 

Evaluating the angular integral of ufg with use of the complete elliptic in- 
tegral of second kind and limiting ourselves to the lowest power in r/R, we 
obtain 
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/(r) = d<t>uUr,4>) ^ ^ (70) 



Introducing this result into Eq. (69), the bending energy becomes 



oo 



i / r/(r)roo(r)dr = ^ / T^^{r)r'dr (71) 







The bending energy leads to a modification of the transverse stability con- 
dition (35). As a consequence, the equation (45) is also modified leading to a 
smaller value for the string width. Applying the spatial rescaling (29) to the 
energy e + Cfe, we obtain 

f'OC f'OC 

/ {dacffrdr = k"^ {di$xd2$frdr 
Jo Jo 

where the second integral on the rhs represents the modification of the con- 
dition (35) caused by bending. We have neglected the contribution to this 
integral due to the term (d^ij))'^. This contribution is logarithmically divergent 
and a cutoff radius of order R must be imposed. Then it turns out that its 
magnitude is less than 1/10 of that due to the second integral on the rhs of Eq. 
(72). Evaluating the integrals in Eq. (72) for the n = 2 Ansatz (41), and using 
the condition (63) for = 1, we obtain the string width in the presence of a 
bend as a function of the loop radius 

If we compare this result with Eq. (45), we see that A;, < A. This is expected 
as the bending energy is decreased on decreasing the width of the string. 

Now, we consider the effect of the bending energy on the stability of the 
loop. Substituting into the rhs of Eq. (71) the quantity Too, and using the 
transverse stability condition (72), we obtain the total energy per unit length 



■ r 

^tot = e + eb= J 



rdr 



+ {d2<i>f + 2{ds(Pf + «2(5,3 + 

+{d2$f + k\Si3 + 523)] + [{ds$f + K^Si2\ jr^^r 



(74) 



Evaluating the integrals in this equation for n = 2 and N = 1, we are led to 
the following result for the total energy per unit length 
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, / \2 , ^ 2, 2 2\2 , 4 2\ 



To simplify numerical calculations, we neglect the small last term on the rhs 
of this eqation and obtain the loop energy in the form 



eioop ^ 2nRv^ An + ? ( 1 + ^) + ^ (76) 



The minimum of this expression occurs for the preferred radius Re satisfying 
the following equation 



;-!7r/v 

+ 



A numerical solution of this equation yields i?g ~ k^. Using this value in Eq. 
(73), we obtain ~ I.Sk^. From these results, we obtain the ratio Re/Xb — 0.81 
showing a substantial enhancement in comparison with Re/X — 0.45 obtained 
from the calculation which disregards the effect of the bending energy. 

Let us compare eioop obtained from Eq. (76) with the numerical results 
obtained in Ref. [17] for the minimum energy configurations characterized by 
Hopf charge Q. According to Ref. [12]. the knot quantum number Q = mn. 
Since 2'jT/k is one period in the x'^ coordinate, we have k = 1/R which implies 
m = N. Thus, the loop with N = l,n = 2 has Hopf charge Q = 2. 

The numerical result for the Q = 2 soliton energy is 835 e.u. [17]. The 
energy unit e.u. = F^/Ae = K,rf /A. By fitting the parameters and e to 
the baryon masses, Adkins, Nappi, and Witten [18] find = 129 MeV and 
e = 5.45. Thus the energy unit used in Ref. [17] is e.u. = 5.9 MeV and the 
Q = 2 soliton energy is E = 4.92 GeV. This should be compared with our result 
for the loop energy 29bKTf ~ 7 GeV. This is consistent with the numerical 
results of Gladikowski and Hellmund [19] who also found that the energy of the 
configuration with A'' = 2, n = 1 is lower than that with A = 1, n = 2. 



9 Gravitational Field of a Circular Loop 

Starting with the Einstein equation in three spatial dimensions, we introduce 
the gravitational potential 0g = /ioo/2 and obtain from the 00-component of 
Eq. (6) a Newton-Poisson equation with the solution 

4>g{S) = -2G [ T^^J^x' (78) 

J \X X \ 
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where t{x') is given in Eq. (13). Developing l/\x — x'\ in a Taylor series, we 
obtain 

(t)g{x) = (t)o{x) + (t)2{x) + - (79) 

where 

0o(x) = -— / T{x')d^x' = (80) 
X J X 

where x is the distance from the loop center. 

The quantity Ma is the total active gravitational mass contained in the 
integration volume. Owing to the rotational symmetry of the loop, there are 
no odd terms in the expansion (79). In what follows, we assume that x ^ R. 
Then (po{x) dominates this expansion. 

For a large loop of radius R, we have Ma = 27: Rnia, where rria is the active 
gravitational mass density per unit length of the string discussed in Sec. 3. 
Using Ansatz (41) on the rhs of Eq. (25), we obtain 

ma = Wi^\j^ + ^-) (81) 

Using in Eq. (81) the equilibrium values Ug — \/5,Re = Nk/VS, and the 
quantity given in Eq.(45) we obtain uia — 497?^. The active gravitational 
mass of the loop is then equal to Ma = 2nRema — IISNkt]'^. 

For a small loop, we need to calculate nia from Eq. (25) with T(r) that is 
modified by the bending energy. The gravitational mass per unit length of a 
rigid string is given by 

POO POO 

fha':^A'K f{r)rdr = 2-k (Tqo + Tn + T22 + %i)rdr (82) 
Jo Jo 

where Tu equals Tu renormalized by the bending energy cb- Thus, we have 

27r / fbordr = 27r / Toordr + (83) 
Jo Jo 



and 



2tt [ fa^rdr = 2w f T^ardr - (84) 
Jo Jo 

Note that the Ihs of Eq.(84) is the net pressure integrated over the cross 
section of the rigid string. The second term on the rhs is the additional pressure 
due to the stress tensor 0-33 given in Eq.(67). This leads to the identity 

poo pIlT poo 

rdrj^ a33{r,^)d^c,^ Too{rydr = (85) 

Next, we consider the remaining terms of Eq. (82). Using Eqs. (9) and (10), 
we have 
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rdr 



(86) 



277 (fn + T22) rdr = nr]^ j [k^Si2 - {d^^f 

For a flexible string, the rhs of Eq.(86) vanishes owing to the stability con- 
dition (35). For a rigid string, this condition is replaced by Eq.(72) so that the 
rhs of Eq. (86) is nonzero and given by 



27r / {Tii+T22)rdr 
lo 

v2 



nrj 

IR2 



+ (52^ +k'(5i3 + 523) 

Using Eqs. (11), (83), and (84), we have 



r^dr 



(87) 



27r 



foa + f 33] rdr = nr]^ 



{d3^^ + K\Si3 + S23) 



rdr 



Introducing Eqs. (87) and (88) into (82), yields 



(88) 



'4i?2 Jo 



W (530)' + k'(5i3 + 523) 

Jo L 

{di$f +{d2$f + {3^3 + 823) 



rdr 



rdr 



(89) 



Performing the integrations over r with use of Ansatze (2) and (41), Eq. 
(89) yields 



(- 



+ ■ 



+ 



2^2 



(90) 



From the numerical solution of Eq. (77), we have fc^A^ = Xl/Rl = 1.5, /c^k^ ~ 
1, and Aj/k^ ~ 1.5. Introducing these values into Eq. (90), we have fha — 34?7^. 
Consequently, the active gravitational mass of the iV = l,n = 2 loop. Ma = 
2'jTRefha, for the nucleon scale is given by Ma — 214k?7^ ~ 5 GeV. This is to be 
compared with the net mass of the = l.n = 2 loop which according to Eq. 
(76) amounts to cioop — 295k?7^ ~ 7 GeV. The decrease of Ma relative to eioop 
is to be attributed to the rigidity of the small loops. 



10 Knots and Dark Matter 

Knowing that the twisted Baby Skyrmion string and loops produce a nonzero 
long range gravitational acceleration, prompts us to look into the cosmological 
implications of the model studied in this paper. In what follows, we focus on 
the question if the smallest rings (knots) could be a candidate for the cold dark 
matter. 
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Spcrgcl and Steinhardt [3] proposed a scenario in which the dark matter 
particles are self-interacting so that they have large elastic scattering cross- 
section but negligible annihilation with a lifetime larger than the present age 
of the universe. Wc note that the model described by Lagrangian (1) is also 
self -interacting owing to the presence of the Skyrme term that is fourth order 
in the field variables. As pointed out in Ref. [12], this term is responsible for 
the topological stability of the rings coming from linking the two color magnetic 
fluxes. In the scenario of Ref. [3], the mean- free path of the particles should 
be larger than 1 kpc but smaller than about 1 Mpc. This implies a range 
8 X 10~(^^~^^) cm^/GcV for the ratio a jm. where a is the scattering cross section 
and m is the mass of the dark matter particles. 

We now examine this ratio for the nucleon scale knot with A'' = l,n = 1. 
According to Ref. [17], the numerical mass for this knot is m ~ 504 e.u. ~ 3 
GeV. Then the corresponding Compton wavelength is Ac = 4 x 10"^'* cm. For 
a rough estimate of the scattering cross-section, we follow Ref. [3] and write 
(T = 47ra^ where the scattering length a ~ lOOAc.Then cr ~ 2 x 10~^^ cm^, 
and cr/m ~ 6.6 x 10~^^ cm^/GeV, a value that is within the range, ajm ~ 
8 X io-(25-22) cmVGeV prescribed by Ref. [3]. 

Next, we consider the electroweak knot for which the energy \init is much 
larger than that for the nucleon scale (e.u.~ 5.9 MeV). Using the relations 
r\'^^? = Ijg^ and My^ ~ •qgj'l where g = ZTicj^ is the coupling constant, and 
Mw ~ 80 GeV is the mass of the W boson, we obtain rf'i^ ~ 5.5?ic and n ~ 
1.74 X 10~^^cm. Thus, the electroweak energy unit is (e.u.)eu; — rfn/A ~ 155 
GeV. Consequently, the mass of the electroweak knot, with N = l,n = 1[17], is 
m ~ 504(c.u.)eu' — 7.8 x 10^ GeV. Taking a =■ A.-kc? where a ~ IO^k, we obtain 
ajm ~ 5 X 10^^2 cm^/ GeV, a value that is about 10"'' times the lower limit 
of the prescribed range [3]. Thus, owing to its large mass, the A' = l,n = 1 
electroweak knot cannot be a candidate for the dark matter particle in the 
scenario of Ref. [3] . 

The circular loops studied in Sees. 7 and 8 are classically stable implying 
that the topological invariant A" is conserved. However, this is not strictly true 
quantum mechanically, since there is a possibility of A^ — >■ A^ — 1 transitions due 
to quantum tunneling. This process involves migration across the loop of the 
flux line that is associated with the current along the loop. We now consider 
the circular A^ = 1, n = 1 loop decaying into the A^ = 0, n = 1 loop. To estimate 
the decay rate, we use the WKB formula 1/t ~ wo exp {—Se) where coo is the 
rate of hitting the barrier and 5'_e — (2K/h)\/2MAe. Here M is the mass of 
the migrating line and Ae ~ etot^ is the energy of a loop segment of width k. 
For A' = l,n = 1, we obtain uiq ^ c/{s/ttk), M ~ nrfn/c? and Ae ~ 2QrfK, 
implying Se — 52r]'^ k"^ / (he) . 

For the nucleon scale loop, we have rj'^K^ ~ 7.1 x lO^^^^^ yielding Se — 3.7. 
With K ~ 4 X lO-^'^cm, the decay time is t ~ 2 x 10"^^ exp (Se) s ~ 8 x lO'^^s 
which is extremely short in comparison with the age of the universe 10"'^''s. This 
excludes the A^ = l,n = 1 knot from being the dark matter candidate within 
the scenario of Ref. [3]. 

This prompts us to examine the stability of the N = loops. Since the anal- 
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ysis of the iV = 0,n = 1 case is complicated by the presence of logarithmically 
divergent integrals, we confine ourselves to the = 0, n = 2 case. 

First, we address the transverse stability condition for this case. Letting 
m = A'' = in Eq. (72), and evahiating the integrals with the Ansatz (41), we 
obtain ~ 2Rk. Next, we consider the loop energy ejoop — '^'^R^tot, where 
etot is given by the expression (75) with m = and ~ 2Rk. In this way, we 



There is a remarkable similarity of this result to the expression (61). We see 
that the third term in the parenthesis of (91) now plays the role of the centrifugal 
energy which stabilizes the loop against a radial collapse. The expression (91) 
has a minimum at the preferred radius given by Re — Using this result 

in Eq. (91), we obtain eioop — ISOt^^k. For the nucleon scale loop, this implies a 
mass ~ 3 GeV. The A'' = 0, n = 1 loop is expected to have a lower mass owing 
to the decreased kinetic energy of the currents (see Eq. (19)). With a mass that 
is smaller than 3GeV, the ratio a /mis expected to be within the range required 
by the dark matter scenario of Ref . [3] . 

We should now inquire about the gravitational mass of the A'' = 0,n = 2 
loop. Letting A^ = in Eq. (90), we have m ~ TT^jy^Af Using A^ ~ 2ReK 
and i?g ~ K^/3, we obtain from this result fria — Mrj"^. Then the gravitational 
mass of the N = 0,n = 2 loop becomes Ma = 2'KRe'fha — 12AffK. Hence, for the 
nucleon scale, rfn ~ 23.6 MeV, we have Ma — 2.9 GeV. Thus the gravitational 
mass of the N = Q,n = 2 loop is practically equal to the entire mass of the loop 
eioop — 3.1 GeV. Consequently, these loops could play a role of the dark matter 
particles in the scenario of Ref. [3] . 

Although, the electroweak knots are too heavy yielding a ratio a/m that is 
seven orders magnitude below the range given in Ref. [3], they may still be of 
cosmological relevance. We recall that in the electroweak scale, the parameters 
ifn^ ~ 5.5?ic and k ~ 1.7 x 10~^^ cm so that the action Se — 286. Thus, the 
decay time of the N = l,n = 1 electroweak knot is r ~ 10^^^ s which is much 
longer than the age of the universe. This is evidently due to the high energy 
cost for the intersection of the flux lines involved in the N ^ N — 1 transition. 

For the same reason, the long electroweak twisted Baby Skyrmion strings are 
prevented from intcrcommuting. This could lead to the formation of a frustrated 
string network. In an interesting paper, Bucher and Spergel [20] proposed that 
a solid dark matter component could arise from frustrated networks of non- 
Abclian cosmic strings. Owing to the finite shear modulus, this solid component 
can have a negative pressure without instabilities that are expected for a perfect 
fluid. Perhaps, the non-Abelian strings studied in the present paper could play 
a role in the theory of Ref. [20] . 



obtain 




(91) 
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